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Abstract. The influence of intrinsic channel noise on the spontaneous spiking 
activity of poisoned excitable membrane patches is studied by use of a stochastic 
generalization of the Hodgkin-Huxley model. Internal noise stemming from the 
stochastic dynamics of individual ion channels is known to affect the collective 
properties of the whole ion channel cluster. For example, there exists an optimal size 
of the membrane patch for which the internal noise alone causes a regular spontaneous 
generation of action potentials. In addition to varying the size of ion channel clusters, 
living organisms may adapt the densities of ion channels in order to optimally regulate 
the spontaneous spiking activity. The influence of channel block on the excitability of 
a membrane patch of certain size is twofold: First, a variation of ion channel densities 
primarily yields a change of the conductance level. Second, a down-regulation of 
working ion channels always increases the channel noise. While the former effect 
dominates in the case of sodium channel block resulting in a reduced spiking activity, 
the latter enhances the generation of spontaneous action potentials in the case of a 
tailored potassium channel blocking. Moreover, by blocking some portion of either 
potassium or sodium ion channels, it is possible to either increase or to decrease the 
regularity of the spike train. 

PACS numbers: 05.40.-a, 87.10.+e, 87.16.-b 

Submitted to: Physical Biology 
1. Introduction 

Following the study of Hodgkin and Huxley pQ , most of the models of axons have treated 
the generation and propagation of action potentials using deterministic differential 
equations. Since the work of Lecar and Nossal [2] it became clear, however, that 
not only the synaptic noise but also the randomness of the ion channel gating itself 
may cause threshold fluctuations in neurons j3J. Therefore, channel noise which stems 
from the stochastic nature of the ion channel dynamics must be taken into account 
jsSj. It impacts such features as the threshold to spiking and the spiking rate itself 
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jH 03 E3 13 El El HO], the anomalous noise-assisted enhancement of transduction of 
external signals, i.e. the phenomenon of Stochastic Resonance [HI El El El and 
the efficiency for synchronization [15J, to name but a few such phenomena. The origin 
of the channel noise [3] is basically due to fluctuations of the mean number of open 
ion channels around the corresponding mean values. Therefore, the strength of the 
channel noise is mainly determined by the number of ion channels participating in the 
generation of action potentials. Interestingly, there is an optimal patch sizes for which 
the spike production becomes more regular fTJ E]- The objective of this work is to 
investigate how the regularity of the spiking can possibly be controlled for a given 
membrane patch size. Toxins like tetraethylammonium (TEA) and tetrodotoxin (TTX) 
allow to reduce the number of working potassium or sodium ion channels, respectively, 
for an extended period of time ^Hj- Moreover, the densities of ion channels can be 
adapted by the living cell also dynamically (TBI Ej over an extended time span which 
is needed, e.g., to express the required ion channel proteins in the membrane ^7j. The 
effect of blocking the ion channels entails several different tendencies at the same time. 
For sodium ion channels, the reduction of the density of channels results in an increase 
of the activation barrier towards excitation from the resting state and, therefore, in 
the reduction of neuronal activity. On the other hand, however, the corresponding 
channel noise component will also be increased due to the reduction of the absolute 
number of ion channels in the membrane patch. As a consequence, the increased channel 
noise will help to overcome the activation barrier and to initiate spontaneous spikes. A 
reduction of the density of potassium channels will on the contrary generally result in a 
lowering of the activation barrier (an increase of the excitability - see, e.g., in ^1) and, 
simultaneously, in an increase of the recovery time which should favor longer interspike 
time intervals. Moreover, the reduction of the total number of potassium ion channels 
will also increase the corresponding channel noise component which is expected to lead 
to an increased variability of the refractory period. Which of these various concurrent 
effects will dominate is not clear a priori. This depends both on the ion channel densities 
and on the size of the studied membrane patch. It is the main objective of this paper 
to study this highly nontrivial, subtle issue in a stochastic model which extends the 
Hodgkin-Huxley model of neuronal excitability. 

2. A stochastic Hodgkin-Huxley model 

According to the Hodgkin-Huxley model the dynamics of the membrane potential V , 
measured throughout this work in mV is given by: 

C^V + G K (n) {V - E K ) + G Na (m, h) (V - £ Na ) + G L (V - E L ) = . (1) 

In Eq. (0), C = 1/iF/cm 2 is the capacity of the cell membrane. Furthermore, 
-^Na = 50 mV, i?K = — 77 mV and E-^ = — 54.4mV are the reversal potentials 
for the potassium, sodium and leakage currents, correspondingly. While the leakage 
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conductance is assumed to be constant, G L = 0.3mS/cm 2 , the potassium and sodium 
conductances read: 

Gk(u) = <?r x *K n\ G Na (m, h) = x^ m 3 h , (2) 

where g^ ax = 36mS/cm 2 and g™ x = 120mS/cm 2 denote the maximal conductances 
(when all the channels are open). In Eq. (j2J) we introduce the factors xk and x^ a 
which are the fractions of working, i.e., non-blocked ion channels, to the overall number 
of potassium, N K , or sodium, iV Na , ion channels, correspondingly. These factors are 
confined to the unit interval. Experimentally, they can be controlled by adding cell 
toxins like tetraethylammonium (TEA) and/or tetrodotoxin (TTX) which completely 
block and disable potassium or sodium ion channels, respectively [18 . 

While the gating variables n, m and h describes the mean ratios of the open gates 
of the working channels, the factors n 4 and m 3 h are the mean portions of the open ion 
channels within a membrane patch. This follows from the fact that the gating dynamics 
of each ion channel is assumed to be governed by four independent gates each of which 
can switch between an open and a closed conformation. The voltage-dependent opening 
and closing rates a x (V) and /3 X (V) (x = m, h,n), read [TT| ITS]: 
n ,s 0-1 (V^ + 40) 

amiV) = l-exp[-(V + 40)/10] ' (3fl) 
(3 m (V) =4 exp[-(V + 65)/18], (36) 
a h (V) = 0.07 exp[-(V + 65)/20], (3c) 

P h (V) = {l + exp[-(V + 35)/10]}-\ (3d) 

0.01(^ + 55) . . 

aniV) = l-exp[-(V + 55)/10]' (3e) 

p n (V) =0.125 exp[-(V + 65)/80]. (3f) 

Fox and Lu [B] have extended the Ho dgkin- Huxley model by taking into account the 
fluctuations of the numbers of open ion channels around the corresponding mean values. 
Within a corresponding stochastic description, the gating variables become stochastic 
quantities obeying the following Langevin equations: 

^■x = a x (V) (1 - x) - p x (V) x + £ x (t), x = m,h,n, (4) 

with independent Gaussian white noise sources £ x (t) of vanishing mean. For an excitable 
membrane patch with N^a, sodium and iVx potassium ion channels the noise correlations 
assume the following form: 

<u*)uo> = r a 7^ ] fi V (L ^ " ^ ■ ^ 

N Nsb i Na [a m {V) + f3 m (V)\ 

iVNa a; Na [oth{V) + p h (V)\ 

N K x K [a n {V) + p n {V)\ 
The overall numbers of involved potassium and sodium ion channels are re-scaled by 
iNa and xk, respectively, in order to disregard the blocked channels which do not 
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Figure 1. The bifurcation diagrams for the undriven original Hodgkin-Huxley model 
are plotted versus the ratio of intact sodium channels XNa, Fig- (a), and the ratio 
of intact potassium channels xk, Fig. (b), respectively, as bifurcation parameters. 
The solid line in Fig. (a) denotes the stable fixed point solution. For potassium 
channel blocking, Fig. (b), the stable fixed point solution (single solid line) becomes 
unstable with decreasing Xk, and a stable spiking and oscillatory solution arises (the 
two solid lines refer to the maximal and minimal membrane potential). The dashed 
lines correspond to instable solutions. 



contribute to the channel noise. With the assumption of homogeneous ion channel 
densities, pNa — 60 /im~ 2 and px = 18 pm -2 , the ion channel numbers are given by: 
-^Na = PNaS*, Nk = PkS , with 5" being the size of the membrane patch. The number 
of working ion channels, i.e. the size of the excitable membrane patch S, respectively, 
determines the strength of the fluctuations and thus the channel noise level. With 
decreasing patch size, i.e. decreasing number of ion channels, the noise level caused by 
fluctuations of the number of open ion channels increases, cf. Eq. (I5aj) - (p)cj) . It is worth 
noting that the Ito-Stratonovich dilemma of the interpretation of the studied system 
of stochastic differential equations does not appear since each of the noise sources £, x (t) 
does not depend explicitly on the state of the corresponding variable x ffiH I2T]. 



3. Poisoning in the deterministic Hodgkin-Huxley model 

Before discussing the impact of channel noise on the spontaneous spiking activity 
of poisoned membrane patches we first consider the role of channel toxins on the 
excitability and resting potential of the original Hodgkin-Huxley model which neglects 
the fluctuations of the number of open ion channels. Then, the equations for the gating 
dynamics [T] reads: 

■^-x = a x (V) (1 - x) — (3 X (V) x, x = m,h,n, (6) 

with the opening and closing rates given by Eq. (|3ajl-(|3j|). Hence Eqs. (JT]), (J2J) and 
(JHJ) form a deterministic Hodgkin-Huxley model which takes poisoning into account. 
Remarkably, this set of equations corresponds to the limit of infinitely large numbers of 
ion channels within the stochastic generalization which was introduced in Sect. |21 

Fig. ^depicts the bifurcation scenario which is derived from the deterministic model 
for the case of poisoning the membrane cluster. While poisoning the sodium channels 
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causes only a small, practically negligible variation of the resting voltage, see Fig. C^a). 
Obviously, the poisoning of sodium channels can - upon neglecting the role of channel 
noise - not induce spiking events. In clear contrast, a reduction of the number of 
working potassium channels changes dramatically the qualitative behavior of the spiking 
activity, see Fig. QJb). Upon decreasing the number of intact potassium channels a sub- 
critical Hopf-bifurcation takes place - a stable spiking and oscillatory solution arises and 
the stable, non-spiking solution becomes unstable. Moreover, there is region between 
0.549 < ik < 0.636, where a stable spiking and a stable fixed point solution coexists. 
With further reduction of the ratio xk the oscillatory spiking solution loses stability, a 
sub-critical Hopf-bifurcation takes place and a stable fixed point solution arises again. 
Once more, a region of bistability within 0.0859 < < 0.1068 is identified. 

4. The mean interspike interval 

The numerical integration of the stochastic generalized Hodgkin-Huxley model, cf. 
Eqs. ffT|- (j5cjl . is carried out by the standard Euler algorithm with a step size of 1 /j,s. 
The Gaussian random numbers are generated by the "Numerical Recipes" routine 
ran2 using the Box-Muller algorithm [22] • To ensure the confinement of the gating 
variables between (all gates are closed) and 1 (all gates are open) we have implemented 
numerically reflecting boundaries at and 1. The occurrences of action potentials are 
determined by upward crossings of the membrane potential V of a certain detection 
threshold. Due to the very steep increase of membrane potential at firing the actual 
choice of the detection threshold does not affect the results. In our simulations the 
spontaneous spikes are found by upward crossings at zero threshold voltage. The 
occurrences of action spikes t iy i = 1, N form a point process. 

While in the original, deterministic Hodgkin-Huxley model the action potentials 
occur only for a certain external current stimulus, the intrinsic channel noise initiates 
spontaneous spikes (3 IH El El El E! • The mean interspike interval, i.e., 



with t — 0, becomes a function of the patch size S. The inverse mean interspike 
interval defines the spiking rate. With increasing noise level or decreasing patch sizes S, 
respectively, the spike production increases and thus the mean interspike interval (T) 
decreases and can approach the refractory time [TT1 112j . 

In Fig. El we depict the mean interspike interval (T) against the fractions x^ a and 
Ik of working sodium channels or potassium ion channels, respectively, for different 
patch sizes S. A reduction of the number of working sodium ion channels, i.e. smaller 
^Na-values increases the strength of channel noise which is caused by the stochastic 
behavior of the sodium ion channels. Because the channel noise, which is induced by the 
sodium channels, is seemingly mainly responsible for the initiation of action potentials 
from the rest potential one might expect that a reduction of the number of sodium 
ion channels could then lead to more spikes. Concurrently, however, a reduction of 
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Figure 2. Mean interspike interval for poisoning: The dependence of the mean 
time between two subsequent, spontaneous action potentials versus (a) the ratio of 
intact sodium channels and (b) versus the fraction of active potassium ion channels is 
shown for four different patch sizes: S — 1 ^m 2 (solid line), S = 4/xm 2 (dashed line), 
S = 16 /im 2 (dotted line), and S = 64 /im 2 (dashcd-dotted line). 



the number of working sodium channels causes a diminishment of the maximal sodium 
conductance. Given this competition between these two mechanisms it is the latter effect 
that dominates and, consequently, causes an increase of the mean interspike interval, cf. 
Fig. El (a). 

The reduction of potassium conductance by poisoning the potassium channels 
changes dramatically the dynamics of the original, deterministic Hodgkin-Huxley model 
which ignores the impact of channel noise. With decreasing potassium conductance 
a stable spiking solution emerges at xk = 0.636. The stable silent solution coexists 
with the stable spiking solution until xk = 0.549. Depending on the initial value the 
deterministic system spends the time either in the silent state, or on the limit cycle for 
0.549 < xk < 0.636. Consequently, the noisy system undergoes noise-induced stochastic 
temporal transitions between both stable manifolds. The motion on the spiking limit 
cycle becomes randomized. For xk < 0.549, the stable fixed point disappears and the 
dynamics takes place on a stochastic limit cycle manifold. Then, the noise does not 
seem to affect strongly the mean period of the stochastic cycling. For different patch 
sizes this mean period (T) is about the same, cf. Fig. EJb). 

However, before the Hopf-bifurcation occurs, the increase of the potassium channel 
noise can enhance the mean frequency of the stochastic spikes occurrences. Here, the 
mean interspike interval (T) can be drastically decreased with a reduction of working 
potassium channels xk, cf. Fig. |2^b). This behavior is rather different from that known 
before. It is generally appreciated [3J |Hj that the fluctuational increase of the number 
of open sodium channels is the reason for the spontaneous spike generation. The role of 
the potassium channel noise is usually thought to be negligible. However, the decrease 
of the potassium conductance is known to enhance the excitability and it can even 
induce spikes deterministically (see before). Therefore, the fluctuational decrease of the 
number of open potassium channel can also induce spikes. Note in this respect that 
about one third of the potassium channels is open at the resting potential. This is 




Figure 3. Same as in Fig. [21 but for the coefficient of variation CV, cf. Eq. JHJ. By 
addition of TEA or TTX the regularity of the spiky output can both be increased and 
decreased. 



contrary to the sodium channels which are almost all closed at the resting potential. 
The described alternative mechanism of the fluctuational spike generation can be clearly 
seen at work in Fig. Ef b) where the sodium channel noise remains unchanged for each 
plotted curve. The effect strongly depends, however, on the size of the membrane patch. 
For a sufficiently large membrane patch (top curve in Fig. 2(b)) the discussed effect is 
indeed very strong. However, when the patch size becomes small (bottom curve in Fig. 
2(b)) the poisoning of potassium channels does not cause much effect since the stochastic 
dynamics is dominated in this case by the channel noise component stemming from the 
stochastic dynamics of sodium channels. These qualitative features are displayed in 
Fig. lib). 



5. Controlling the Coherence of Poisoned Spiking Activity 

We next address the regularity of spontaneous action potentials. A proper measure is 
the coefficient of variation, CV, a measure of coherence, which is given as the ratio of 
standard deviation to the mean value: 



CV-J^-, (8, 

where (T 2 ) := -rrJ2(ti — ti-i) 2 is the mean-squared interspike interval. For a fully 
disordered point process (the case of Poisson process) the coefficient of variation CV 
assumes the value CV = 1, while for more ordered processes it assumes smaller values 
and for a deterministic signal it vanishes. In previous studies it has been demonstrated 
that CV exhibits a distinct minimum for an optimal patch size S ~ 1 /mi 2 at which the 
spiking is mostly regular ^] . This phenomenon has been termed intrinsic coherence 
resonance. 

In Fig. the coefficient of variation is plotted vs. the fractions iNa (a) and i K 
(b). Any addition of TTX which blocks sodium ion channels leads to an increase of the 
CV, cf. Fig. Efa). In contrast, toxins which disable potassium channels yield a rise of 
the regularity, even beyond the level which can be reached for an optimal patch size 
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with the unmodified density of the ion channels, cf. Ref. jTTj. This result is due to the 
fact that the lowering of the potassium conductance level promotes a stable oscillatory 
spiking solution in the deterministic model. The channel noise, on the contrary, tends 
to disturb the regular spiking. Nevertheless, the former effect dominates resulting in the 
observed behavior for the CV. In the case of optimal patch size, where the phenomenon 
of intrinsic coherence resonance occurs, the poisoning of potassium channels does not 
produce a strong effect. This is in a clear contrast with the poisoning of sodium channels 
which destroys the coherence resonance (bottom curve in Fig. 3(a)). 

6. Conclusion and outlook 

Our study of the stochastically generalized Ho dgkin- Huxley model reveals the possibility 
to manipulate the response of a spiking membrane patch by adding toxins which 
selectively block ion channels. For example, by a fine-tuned addition of tetrodotoxin 
(TTX) a certain portion of sodium ion channels could be experimentally disabled on 
purpose. This in turn results in a reduction of the spontaneous action potentials and 
causes a more irregular production of spikes. On the other hand, the addition of 
tetraethylammonium (TEA) can be used in order to block potassium ion channels. 
This causes a surprising increase of the spiking activity (i.e. a decrease of the mean 
interspike interval) and yields in turn a more regular spontaneous spiking coherence. 
These characteristic features are expected to impact as well the behavior of biological 
"Stochastic Resonance" , i.e. the phenomenon that the application of an appropriate dose 
of noise can boost signal transduction [TTT^ fT^] and, also, the phenomenon of "Coherence 
Resonance" [22] m oscillatory or excitable biological entities [TT] . 

We share the confident belief that our study of the tailored control of channel 
noise via channel blocking in an archetypal model of excitable biological membranes 
provides some insight into the underlying principles and mechanisms and thus will 
motivate further studies of more realistic models of real neurons where such channel 
noise phenomena do play an essential and constructive role. 

Acknowledgments. This work has been supported by the Deutsche Forschungs- 
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Glossary 

Action potential / spike. A rapid change of the trans-membrane electrical potential in 
the excitable cell membrane. 

Bifurcation. A qualitatively change in the topology of the attractor-basin phase portrait 
under a small variation of a parameter within a nonlinear system. 

Channel noise. The fluctuations of the number of open ion channels. 

Coherence Resonance. Noise-induced improvement of the regularity of the system's 
output. 



Effect of channel block on the spiking activity 



9 



Excitable membrane. Cell membranes that are capable of rapid changing their trans- 
membrane electrical potential in form of spikes. 

Gates / gating. Closing and opening function of an ion channel. This usually refers 
to the movement of protein structural elements of the channel that occludes the 
channel pore. 

Gating variable. A variable which accounts for the fraction of open channel gates of a 
certain type. 

Gaussian white noise. Stochastic trajectories whose distributions are uncorrelated in 
time and are Gaussian-distributed. 

H op f -bifurcation. A bifurcation of a fixed point to an oscillatory solution. When 
the oscillatory solution is unstable and exists at subcritical values of the control 
parameter, we term it a subcritical Hopf-bifurcation. Conversely, when the 
oscillatory solution is stable and exists at supercritical values of the control 
parameter, a supercritical Hopf-bifurcation is observable. 

Ion channel. A protein folded in the cell membrane which enables the transport of 
specific ions through the membrane. 

Ito-Stratonovich dilemma. Interpretation problem which arises in the context of 
Langevin equations in case of multiplicative Gaussian white noise. 

Langevin equation. An equation of motion that describes the temporal evolution of a 
variable which is subjected to noise acting on the system. 

Limit cycle. An attracting set of a nonlinear system to which trajectories converge. 
Upon the limit cycle these trajectories are closed and periodic. 

Manifold. A topological space wherein a stable solution of a nonlinear system remains. 

Membrane patch. A portion of the cell membrane with a certain size. 

Membrane potential. The trans-membrane electrical potential. 

Opening / closing rate. Under the two-state assumption for the gate dynamics (open 
and closed) these give the probability per time unit for transitions between the two 
states of a single gate. 

Poisson process. A stochastic process which is memory-less (Markovian) with 
exponentially distributed waiting times between two successive events. 

Recovery time. The time needed for relaxation to the resting state after an excitation 
occurred. Under physiological conditions an initiation of an action potential within 
this time span is not possible. 

Reversal potential. The electrical potential for which the trans-membrane flux of specific 
ions vanishes. 

Rest potential / rest state. The equilibrium position of the trans-membrane potential of 
a certain membrane patch. 

Spontaneous spiking activity. The occurrence of action potentials which are not initiated 
by an externally applied stimulus. 
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Stochastic Resonance. A anomalous, noise-assisted enhancement of transduction of 
weak (deterministic or stochastic) signals. 
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